arXiv:chao-dyn/9808020vl 23 Aug 1998 


On the Heat Transfer in Rayleigh-Benard systems 


R. Benzi 1 , F. Toschi 2,3,4 , R. Tripiccione 4 
1 Dipartimento di Fisica, Universita di Tor Vergata 
Via della Ricerca Scientifica 1, 1-00133 Roma, Italy 
2 Istituto Nazionale di Fisica della Materia, unita di Tor Vergata 
3 Dipartimento di Fisica, Universita di Pisa, 

Piazza Torricelli 2, 1-56126, Pisa, Italy 
4 Istituto Nazionale di Fisica Nucleare, sezione di Pisa 
(February 5, 2008) 

In this paper we discuss some theoretical aspects concerning the scaling laws of the Nusselt 
number versus the Rayleigh number in a Rayleigh-Benard cell. We present a new set of numerical 
simulations and compare our findings against the predictions of existing models. We then propose 
a new theory which relies on the hypothesis of Bolgiano scaling. Our approach generalizes the one 
proposed by Kadanoff, Libchaber and coworkers and solves some of the inconsistencies raised in the 
recent literature. 


I. INTRODUCTION 

In this paper we discuss the scaling properties of heat transport in Rayleigh-Benard systems. Although this topic 
has received a lot of attention during the last years, still there are controversial interpretations of experimental results. 
In 1989, Libchaber, Kadanoff and co-workers [[jj have shown that a new and unpredicted scaling range is observed in 
Rayleigh-Benard systems at large enough Rayleigh numbers. They also proposed a physical picture of the experimental 
observations, based on the crucial role played by “plumes”, the coherent structures present in thermal convection, 
in heat transport. This picture was somehow questioned by Shraiman and Siggia, who proposed a rather different 
theoretical approach. Both models make relatively ad hoc assumptions on qualitative features of the fluid motion in 
appropriate sub-domains of a Rayleigh-Benard cell. These assumptions, although reasonable on physical ground, are 
not explicitly justified by the still poorly understood dynamics, in a system of rather complex phenomenology. 

In this paper, we try to build a more dynamical picture of the problem, by establishing a link between the scaling 
properties of heat transport and other well controlled dynamical features encountered in convective turbulence. We 
start presenting and discussing a new set of results based on direct numerical simulation which, in our opinion, are 
able to shed more light on the problem. We then generalize the model of Castaing et al. [jp with a simple ansatz, whose 
basic physical meaning is that the dynamics that control Bolgiano scaling in convective turbulence also modulates the 
scaling of the Nusselt number. Our ansatz leads to well defined predictions that we have verified numerically (and 
could be experimentally checked). The obvious objection that Bolgiano scaling (expected to set in at scales larger 
than the Bolgiano length) can hardly play any rolein heat transfer (where scales of the order of the boundary layers 
are obviously important) can be solved by observing, as already pointed out in |Ii|, that quantities usually regarded 
as globally characterizing the flow (such as energy or temperature dissipation, Bolgiano length) can still be locally 
defined in a convective cell as a function of the vertical coordinate. In particular, Bolgiano length decreases sharply 
in regions close to the boundary layer. 

Our paper is organized as follows. In section || we introduce the equation of motion, the dimensionless parameters 
describing the problem and Bolgiano scaling. In section we review the theoretical models proposed by Castaing et 
al. and Shraiman and Siggia emphasizing the physical assumptions underlying the two different approach. In section 
0 we present our numerical simulations and in section [v|, the most important part of the present work, we discuss 
and justify a simple ansatz which generalizes the model proposed in jp . In section 0 we present our conclusions. 


II. THE PROBLEM 

We consider a fluid in a rectangular cell of horizontal size L and vertical size H. The fluid is heated from below 
and cooled from above by contact with two heat reservoirs. The temperature field T(x,y, z;t) satisfies the boundary 
conditions: 


^ AT 


- AT 


TL=o = T — T| z=ff = T+ — 


(1) 


1 


Furthermore the vertical walls are supposed to be adiabatic (i.e. thermal exchange through these walls are supposed 
to be negligible). We also assume that the fluid flow satisfies the Boussinesq equation of motion (see ||): 
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where g is the acceleration of gravity, a the volume thermal expansion coefficient, v and x kinematic and thermal 
diffusivity, 9 = T — T. The Boussinesq equations of motion are a first order approximation of buoyancy effects and 
are supposed to be accurate if thermal gradients are not too strong. Furthermore the incompressibility condition is 
valid if velocities are much smaller than sound speed. 

The relevant dimensionless parameters in equations and are the Rayleigh number Ra , the Prandtl number Pr 
and the aspect ratio F: 
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For large enough values of Ra , the flow becomes turbulent, i.e. chaotic both in time and in space. One of the basic 
issues that we want to discuss in this paper concerns how much heat is transferred from the bottom boundary to the 
top boundary. To this end it is useful to introduce the dimensionless number: 


_ <«,r> - x^ _ <u/0'> 

- X AT ~ AT 

A H 


( 6 ) 


where u/ is the turbulent fluctuation of the vertical velocity, 9' the (turbulent) temperature fluctuation and (...) 
stands for space and time averages. Nu is called the Nusselt number. It measures the amount of heat transported 
by turbulent fluctuations with respect to the heat transport due to molecular motion ( x H T is the heat flux due to 
conduction if the fluid were at rest). 

We can now state our problem in terms of the dimensionless parameters so far introduced: we want to investigate 
the functional relationship, 


Nu = f(Ra,Pr,T) 


( 7 ) 


To our knowledge, no meaningful effect has been found experimentally on turbulent heat transport due to the 
geometric parameter T. Therefore, we shall neglect, hereafter, the dependency on T in |?]. We shall confine most of 
our discussion to the Pr = 1 case. 

If the fluid were at rest the temperature drop inside the cell would be linear. Because of convection the mean 
temperature profile differs from the purely conductive one. In particular one find that almost all the temperature 
drop occurs in the thermal boundary layers. The thickness of the thermal boundary layer, A, can be roughly defined 
as the thickness of the layer over which the temperature drop is nearly equal to AT /2. An important relation 
(experimentally well verified) allows us to connect A to the Nu number: 


Nu ~ 


H 

~X 


We remark that relation @ becomes rigorously true if one defines A as follows: 
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By means of the relation |8| our problem can be rephrased as the understanding of the scaling of A versus Ra. 

A related issue that we want to address is the role of buoyancy forces effects on the statistical properties of turbulent 
fluctuations for Rayleigh-Benarcl convection. Let us introduce the velocity and temperature difference defined as: 


Sn//{r) = [v(x + r) — v(x)] • r 
ST(r) = T(x + r) - T{x) 
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Following Yakhot |li| we can write: 
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where e is the mean rate of energy dissipation, e = j f )T/ • (diVj + djVi) 2 d 3 x, and N is the mean rate of temperature 
dissipation, N = j f y\ ( d{T ) 2 d 3 x. 

Equations [l^ and |13| replace the well known “4/5” Kolmogorov equation for homogeneous and isotropic turbulence. 
These equations have been derived by assuming that small scale turbulent fluctuations in Rayleigh-Benard system 
are, to first approximation, homogeneous and isotropic (see Yakhot [^2|). The second term on the right hand side of 
|l2| and |l3| represents however a non isotropic, thermally driven, contribution. Therefore |l2| and [l3|, although cannot 
be proven rigorously, can be useful as guidelines for describing the difference, if any, between homogeneous isotropic 
turbulence and thermal turbulence. 


Neglecting intermittency effects, from 13 we have the scaling relation: 

5v(r)5T(r) 2 ~ Nr 


(14) 


There are two physically interesting limit regimes in equation |12p either the first or the second term on the right 
hand side dominates. In the latter case we deduce the following balance: 


5v 3 (r) ~ agST(r)5v(r) ■ r 


By using |h| we obtain: 


Sv 2 (r ) ~ agrST(r) ~ agr 


Nr 
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We can equivalently rephrase eqns. |15j- 16 as scaling laws for the velocity and temperature increments: 

Sv(r) ~ (ag) 2 ^ 5 iV 1 / 5 ?- 3 / 5 
ST(r) ~ (ag)~ 1/5 N 2/5 r 1/5 
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The scaling property defined by [I7| - |l8| is referred to as Bolgiano scaling 0- Using Bolgiano scaling, we can evaluate 
the second term of the right hand side of [l2]. We obtain: 

r' 4 (ST(r')Sv(r')) dr' ~ agST(r)5v(r) ■ r ~ ( ag ) 6 ^ 5 N 3 ^ 5 r 9 ^ 5 (19) 

From |l9]we obtain a consistency condition, requiring that the second term of the r.h.s. in [l2]is larger than er. This 
is true for scales r > Lb, where: 



Lb = 


e 5/4 

TV 3 / 4 (ag) 3/2 


( 20 ) 


Ls is called the Bolgiano scale. For r < Lb, the statistical properties of thermal turbulence should be described 
by Kolmogorov theory of turbulence with scaling: 

Sv(r) ~ e 1/3 r 1 / 3 (21) 

5T(r) ~ N 1/2 e~ 1/6 r 1/3 (22) 

We remark that equation ^2] corresponds to the scaling of a passive scalar while equati on pl| is the usual Kolmogorov 
(1941) prediction. Finally let us notice that for r > Lb the first term on the r.h.s. of |l3| is always greater than the 
second term, which represents a consistency condition. 

The dependence of both N and e on Nu and Ra can be exactly derived from the equation of motion. We have (see 
reference || for details): 
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(23) 


(e) = (Nu - 1 )-Ra 

(N) =Nu-Ra 2 (24) 

H 3 ( ag ) 

Beside any rigorous derivations, we can give a physical meaning to |2^, using the following arguments. Consider 
first the mean rate of temperature dissipation. Since almost all the temperature gradient (AT) is across the thermal 
boundary layer (of thickness A), we can estimate: 

(25 > 

On the other hand, by using equation |23|, we obtain: 

5v(L b ) ■ 6T(L b ) — Nu ■ Ra (26) 

which is equivalent to say that the Bolgiano length L B can be interpreted as the characteristic scale of the eddies 
which transport heat in a Rayleigh-Benarcl system. 

Evidences of Bolgiano scaling 0 and 0 have been reported both in 2D Q] and 3D Q numerical simulations of 
Rayleigh-Benard systems (see |j) for a detailed description including the effects of intermittency). 


III. REVIEW OF PROPOSED SCALING THEORIES ON HEAT TRANSPORT 


In this section we review some theoretical models, proposed in the past, to derive the scaling properties of Nu vs. 
Ra (for a review see also jj|). Our emphasis is on the physical assumptions underlying the models, rather than on 
rigorous derivations from the equations of motion. 

We first discuss three different arguments leading to the old (and currently experimentally disproved) scaling 
relation Nu ~ Ra 1 ^ 3 . Such relation was predicted by many authors, among them Malkus |19|2(J, Priestley [ fjT|| , 
Howard [ ^2|j23| , Spiegel |p4[ . 

The first argument assumes that the boundary layer is marginally stable, i.e. that the effective Ra( A) number 
computed at the thermal boundary layer thickness, A, is equal to a critical Rayleigh number Ra c , independent of 
Ra. Since by definition, Ra = agATH 3 /v\ : we get Ra{ A) = Ra ■ (A /H) 3 and hence Ra = Ra c (Hj A) 3 . Therefore, 


A ~ H I 


1/3 


Using 0 we finally obtain the result Nu ~ Ra 1 / 3 . 


Another way to reach the same result is the assumption of the independence of the heat flux from the height of the 


cell, H. Supposing a scaling of the form Nu ~ Ra 7 , the Nusselt number is just 


^ agATH 3 y ^ 


7 is required to 


be 1/3, for the heat flux F be independent of H. From this argument we understand that every model which leads 
to a 1/3 exponent, does, implicitly, assume a decoupling of the top and bottom boundary layers. 

Finally, we want to point out a rather simple argument proposed in | |l2| leading to the same result. Let us assume 
that, close to the thermal boundary layer, the statistical properties of turbulent fluctuations are not affected by 
buoyancy effect. This implies, that we can use the standard dimensional analysis of Kolmogorov theory. Let us also 
assume that the thermal boundary layer thickness, A, equals the viscous boundary layer thickness, 77 , (i.e. A ~ 77 ). 
Assuming K41 scaling for velocity differences, 6v(r) ~ (er) 1 ^ 3 , from the condition Re ~ 1 ~ Sv (v)v anc j f rom the 
scaling e ~ Nu ■ Ra we obtain Nu ~ Ra 1 / 3 . We remark that this argument is self-consistent, because e ~ Ra 4//3 , 
N ~ Ra 7 / 3 so that L B ~ Ra _1 / 12 i which implies g ~ g -1 / 4 ~ Ra~ x / 3 and A ~ 77 < L B . The above argument can 
also be derived by estimating the energy dissipation e(A) in the thermal boundary layer by the Kolmogorov relation 
e ~ w(A) 3 /A, where u( A ) 2 ~ (agATX). By assuming that most of the energy dissipation takes place in the thermal 
boundary layer and using [24| we obtain Nu ~ Ra 1 / 3 . 

We observe that in all cases, for scales of order A and Ra large enough, A < Lg- Thus Bolgiano scaling should not 
be applied near by the boundary layer. 

In an important paper, Castaing et al. jjj reported for the first time clear evidence that 


Nu 


Ra 7 ; 


7 


~ 0.281 ~ 


2 

7 


(27) 


Castaing et al. (1989) have shown that this scaling law is valid for Ra > 10 6 . Further experimental results 
(Ciliberto ^|) have shown that 7 ~ 2/7 even for lower Ra. The results reported in Castaing et al. have motivated 
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many experimental and theoretical efforts aimed at understanding the physical mechanism leading to the (unexpected) 
scaling [27| Here we review two rather different theoretical model proposed in Castaing et al. [l| and Shraiman-Siggia 


The Castaing et al. flj] theory is based upon the assumption that there exist three layers (hereafter referred to as 
A, B and C) characterized by different physical properties: 


A-layer: the thermal boundary layer, near the bottom and top boundaries, of thickness A and temperature differences 
AT/2. In the A-layer, instabilities generate plumes, of typical size A, which are expelled into the B-layer. 

B-layer: a mixing region, of thickness much greater that A and smaller than the size of the cell H. In this layer 
thermal plumes are accelerated due to buoyancy effect. We shall indicate by ST and Sv the the characteristic size of 
temperature and velocity fluctuations, respectively. 

C-layer: the central region of the cell of size comparable with the size of the system. Velocity and temperature 
fluctuations will be indicated by u c and T c respectively. In this layer thermal plumes are advected with almost 
constant velocity. 


The physical picture behind this theory is the following. Thermal plumes are generated in the A-layer, accelerated 
in the B-layer and advected in the C-layer. The Nusselt number can be estimated at the center of the cell (where the 
heat flux is purely convective) as 

Nu ~ u c T c (28) 

To estimate the velocity fluctuations in the center of the cell the only dimensional relation (ignoring thermal 
diffusivity and kinematic viscosity) is the following 

u c ~ (gaT c ■ if ) 1 / 2 (29) 


The basic assumption of the theory is that in the B-layer the characteristic velocity fluctuation are given by the 
balance between the buoyancy effect of the plume and the viscous effect, while temperature fluctuations are equal to 
the temperature fluctuations carried out by the plumes, namely AT. Furthermore, the velocity fluctuations in the 
B-layer equal the velocity fluctuation in the central region. Thus, we obtain: 


u c ~5u~ gaATA2 ST ~ AT (30) 

V 

From relation ||, ^ 8 |, ^9] and M it follows that Nu ~ Ra 2 / 7 . Let us remark that, due to ||o], A can be regarded as the 
cutoff-scale of velocity fluctuations. One of the most important point in the Castaing et al. theory is the role played 
by the thermal plumes, which are well identified coherent structures observed in the chaotic dynamics of thermal 
turbulence. Equation ^o| is based on the assumption that coherent plumes do exist and are observed in the B-layer 
and set the characteristic velocity and temperature fluctuations in Rayleigh-Benard systems. 

A major criticism on the Castaing et al. model concerns the validity of equation |3(] in the mixing layer B. Indeed, 
as already discussed, equation ^ implies a balance between viscous dissipation and buoyancy force. The buoyancy 
force, however, should be relevant only for scales larger than Lb, the Bolgiano scale. The analysis performed in 
section plj indicates that, for scales close to the thermal boundary layer A, Lb is much greater than A, which implies 
that velocity fluctuations cannot be controlled by the overall strength of the buoyancy force. This implies that the 
thickness of the thermal boundary layer cannot be fixed by the balancing proposed in Castaing et al. 

In order to overcome this criticism, a rather different approach has been proposed by Shraiman and Siggia, whose 
theory is based on the relevant dynamical role played by the mean flow observed in Rayleigh-Benard cells. The onset 
of a mean flow is due to plumes rising from the unstable boundary layer. On the other hand, the mean flow generates 
a viscous boundary layer which, in turns, control the thickness of the thermal boundary layer. It is assumed, therefore, 
that the thermal boundary layer is contained in the viscous layer. As we shall see, the most important assumption in 
the theory is that all energy dissipation is constrained inside the viscous boundary layer. 

The starting point of Shraiman and Siggia theory is that, within the thermal boundary layer, there exists a balancing 
between the mean flow advection of horizontal temperature gradient and vertical thermal dissipation, namely 


dT d 2 T 
~ X ~d^ 


(31) 


The velocity profile is supposed to be linear inside the viscous boundary layer (r is a mean shear which has to be 
determined self-consistently) 
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u 


(32) 
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From H and |3^ we can derive a relation between the thermal boundary layer thickness (A) and the unknown 
parameter r. 


A ST ST 

A 3 - X Lt (33) 


The viscous layer thickness ( 77 ) can be estimated using the requirement 


77 

7 • V 


~ 1 


(34) 


Using the exact relation (e) = Nu ■ Ra and under the assumption that the relevant part of the energy is dissipated 
inside the viscous boundary layer, we obtain 


(e) = Nu ■ Ra ~ — 

T Z 


H 


(35) 


Finally from 33, [m] and |35| we obtain Nu ~ Ra 2 / 7 . 


We want to point out that the two theories are based on two quite different physical pictures of the basic mechanisms 
which control heat transport in Rayleigh-Benard system. In the theory by Castaing et al. buoyancy effects and 
their balance with dissipation control the characteristic size of temperature and velocity fluctuations. According to 
Shraiman and Siggia, on the other hand, buoyancy is responsible to maintain the mean flow in the system which 
dynamically controls temperature and velocity fluctuations. 

In the next section we present further experimental and numerical results which will clarify the physics of the 
thermal boundary layer. 

We close this section remarking that the scaling properties discussed so far, are not supposed to be asymptotic, 


as shown by Kraichnan 1111 in the early sixties. For very large Ra numbers an asymptotic regime should emerge as 
can be understood by the following argument. The maximal velocity which can be reached inside the cell can be 
estimated as: 


U M = (agATH) 1/2 (36) 

Since the maximal rate of energy dissipation associated with this velocity is U'l f /H = e, it follows that Nu ~ Ra 1 ! 2 . 
We notice that for such a scaling regime the Bolgiano scale would became Ra independent. Asymptotic prediction of 
1/2 exponent can also be derived as a rigorous upper bound of heat transport (see Doering and Constantin). 

From the experimental point of view, some evidences of transition to the asymptotic regime have been recently 
reported |u|. 


IV. NUMERICAL RESULTS 

An important insight in the physics of thermal convection can be obtained by direct numerical simulation (DNS) of 
Rayleigh-Benard cells. An obvious limitation of DNS is that the available range in Ra is usually much smaller than 
in laboratory experiments. Nevertheless, DNS are extremely useful in understanding the validity of different physical 
assumptions and, occasionally, in checking the predictions on physical quantities which are difficult to measure in a 
real experiments. In this section we discuss a set of new numerical simulations aimed at understanding the correct 
physics of heat transport in Rayleigh-Benard . 

As already pointed out, numerical simulation are confined to rather small Ra number with respect to those available 
in real experiments. In ^j] the 2/7 scaling exponent was observed for Ra larger than 10 6 . The same scaling exponent 
was reported by De Luca et al. |hJ in 2D numerical simulation of Rayleigh-Benard. 

Here we show that the 2/7 scaling exponent is clearly observed, in 3D, even at rather low Ra number. We have 
performed a number of numerical simulations using a LBE (Lattice Boltzmann Equation) [ p~4| code on a parallel 
supercomputer (]l7j for a 3D Rayleigh-Benard cell of aspect ratio 1. For a detailed description of the numerical code 
used we refer the reader to 0 . For the temperature field, we have imposed adiabatic boundary conditions on vertical 
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walls, while top and bottom walls have been kept at fixed temperature (respectively — AT/2 and AT/2). Velocity 
boundary conditions are free-slip on vertical walls and no-slip (zero velocity) on top/bottom walls. All run has been 
done at Pr = 1 and cover a range of about one order of magnitude in Ra around Ra ~ 10 6 . 

We have performed runs at different Ra in order to measure the Nu(Ra) dependence. We found a clear Nu Ra 2 / 7 
scaling, see Figure [1]. 

The results shown in Figure [1] together with those reported by De Luca et al., clearly indicate that the scaling 
exponent does not depend on the dimensionality of the system. This implies that arguments based on three dimensional 
Kolmogorov theory of turbulence must be ruled out. 

As discussed in the previous section, the viscous boundary layer near the top and the bottom of the Rayleigh- 
Benard cell plays an important role in the model suggested by Shraiman and Siggia. In order to clarify this point, we 
have performed a new set of numerical simulations, with free-slip boundary conditions for the velocity on the top and 
bottom walls. The choice of this kind of boundary conditions allowed us to completely remove the velocity boundary 
layers, while keeping a dissipation of kinetic energy inside the bulk of the cell. In Figure [1] we have also reported the 
scaling of Nu versus Ra for the free-slip case. A clear scaling is observed, again, with slope close to 2/7. 

The numerical simulations discussed above indicate that the viscous boundary layer plays a marginal role in de- 
terming the scaling exponent of heat transport. We want to remark that in the theory proposed by Shraiman and 
Siggia one of the basic physical assumptions is that most of the energy dissipation occurs in the viscous boundary 
layer. Certainly this is not the case in the numerical simulation reported in Figure [1] (free slip boundary condition), 
although the scaling exponent of Nu does not change. Our findings, therefore, seem to rule out the approach proposed 
by Shraiman and Siggia, at least for what concerns the assumption on the energy dissipation. 


V. BOLGIANO LENGTH AND NON-HOMOGENEOUS CONVECTIVE CELL 


We have seen in section III that one of the major criticism on the theoretical model proposed by Castaing et al. 
concerns the relevance of buoyancy force in the thermal boundary layer. The criticism is based on the fact that the 
Bolgiano scale Lb is much larger than A for large Ra number. By using the observed scaling A ~ Nu -1 ~ Ra 2 / 7 , 
together with equation |Ts|. p~i| and ^o] we obtain Lb — Ra ~ 3 / 28 . This result show that the ratio Lb /A increases as 
Ra 5 / 28 . From a physical point of view, the Bolgiano scale represents the scale at which energy is injected in the 
system as ’’potential energy”. Buoyancy force converts this energy in kinetic energy. 

The analysis made in section || was appropriate for a homogeneous/isotropic convective cell. However, Rayleigh- 
Benard convective cell is not homogeneous neither isotropic. We can slightly generalize the analysis of section 0by 
assuming that turbulence in the Rayleigh-Benard is “locally” homogeneous and isotropic. It follows that we must 
interpret the various quantities (as for example the energy dissipation e(z), the Bolgiano length Lb(z) and so on) as 
depending locally on z: the distance from the bottom wall. 

Following this approach, we can introduce a local (but plane-averaged) Bolgiano length as the following: 


L b (z) = 


s(z ) 5 / 4 


N(z) 3 / 4 ( ag) 3 ^ 2 


(37) 


Using direct numerical simulations, we can obtain the behavior of e(z) and Lb(z) in the Rayleigh-Benard. In Figure 
[2] we report, the energy dissipation averaged on horizontal planes as a function of 0 . As it is evident, energy is not 
evenly dissipated inside the cell. In Figure [3] the typical behavior of the Bolgiano length, Lb(z), as obtained from 
definition |37], is shown. 

We remark that while the Bolgiano length grows inside the bulk of the cell (in particular in the center of the cell it 
reaches its maximum, nearly equal to the size of the cell itself) it is rel ativ ely small near the top/bottonr boundaries. 

This observation solves the apparent inconsistency raised in section HI where we noticed that A -C Lb' while the 
global Bolgiano length, Lb, can be much greater than the thermal boundary layer thickness, locally, the Bolgiano 
length, Lb{z) has its minimum value around the boundary layer thickness itself. 

The results shown in Figures [2] and [3] suggest that close to the thermal boundary layer the buoyancy force is 
the dominant effect in the system. Thus, velocity and temperature fluctuations should be described by the Bolgiano 
scaling [ll] and [jj]. On the other hand, the thermal boundary layer can be interpreted as the scale at which dissipation 
becomes relevant with respect to buoyancy force. By using and |l8] and balancing the second term (namely the 
forcing) in equation |T^ with the third term, we can introduce the Bolgiano dissipation length tb- After some algebraic 
computation we obtain 
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r B ~ iV-i/8 


(38) 


r\e is equivalent to the Kolmogorov length 77 for turbulent flows without buoyancy force. In a Rayleigh-Benard cell, 
temperature dissipation, N, is confined in the thermal boundary layer. Therefore, equation ^ 8 | can be used locally 
near the thermal boundary layer assuming that N is the global mean rate of temperature dissipation. 

Looking again at the results shown in Figure [2] and [3], we are tempted to assume that in Rayleigh-Benard 
turbulence the thickness of the thermal boundary layer adjusts itself in such a way that it becomes equal to the 
Bolgiano dissipation length, i.e.: 


\ — r B 


(39) 


Using equation |38| and 24 we immediately obtain Nu ~ Ra 2 /'. Our ansatz |39| generalizes the approach proposed by 


Castaing et al. by using the basic equation |l2j and 13. 

By using the definition of L B (z) and the Bolgiano scaling 111 and H we can compute the value of the Bolgiano 
length at scale A. We estimate the rate of energy dissipation as e(A) ~ Sv ^ 
finally obtain: 


Using 24 for the estimate of N, we 


L b { A) ~ 1V - 1 / 4 A -1 

From equation we immediately see that the ansatz |39] implies: 

L b ( A) A 


(40) 


(41) 


Equation |d] should be considered a prediction of the theory so far discussed. We show in Figure [4] the values of 
L b { A) and A as obtained in our numerical simulations. As one can see, the numerical results are quite consistent with 
the prediction |d]. 

Inside the thermal boundary layer, we can assume, following Shraiman and Siggia, that the velocity profile is linear 
in z and that the advection of the mean flow balances the thermal dissipation as expressed in equation [nj By using 
the same approach suggested by Shraiman-Siggia, we obtain that the mean shear inside the thermal boundary layer 
should be proportional to A -3 . This implies that the mean rate of energy dissipation e B L ., integrated inside the 
thermal boundary layer, is proportional to A -5 . Therefore we finally obtain: 

£b.l. ~ R& ^ e TO T (42) 


Equation ^2] is consistent with the numerical simulations, see Figure [5] where we show e B L . plotted against e TO T in a 
log-log plot. Equation [T^ puts bound on the range of Ra where the 2/7 scaling regime is observed. More precisely, a 
transition to the asymptotic scaling predicted by Kraichnan should be observed for values of Ra such that e B L . ~ e TO T- 
According to Figure [5], the critical Ra number for this transition is predicted at Ra ~ 10 n±1 for Prandtl number of 
order 1 . This prediction is in reasonable agreement with recent experimental results found by Chavanne et. al. |p~5|. 


VI. CONCLUSIONS 

The most important result shown in this paper is that the observed scaling of Nu versus Ra in Rayleigh-Benard 
systems can be explained by assuming that the thickness A of the thermal boundary layer is controlled by (and also 
approximately equal to) the Bolgiano dissipation scale r B , i.e. the scale at which buoyancy forces are balanced by 
the dissipative effects. In order to justify our assumption we have introduced a local Bolgiano scale L B (z) based on 
the local values of energy and temperature dissipation. By using direct numerical simulations of Rayleigh-Benard 
turbulence, we have found that L B (z ) is rather small near the thermal boundary layers and becomes equal to the 
cell size near the center of cell. This is a key point in our analysis because it allows to use Bolgiano scaling in order 
to compute the dissipation scale. Previous investigations started with the observation that L B is much larger than 
A and, hence, for z ~ A Bolgiano scaling should be ruled out. Using the assumption that A ~ r B , we predict that 
L b ( A) ~ A. Such a relation has been verified numerically. Our findings support and, in some sense, generalize the 
model proposed few years ago by the Chicago group Q. Let us also remark that our results agree quite well with the 
observed Bolgiano scaling in direct numerical simulations of Rayleigh-Benard systems. 

We have not discussed the dependence of Nu on the Prandtl number Pr. Both the model proposed by Castaing et 
al. and Shraiman-Siggia predicts a Pr dependency which scales as Pr~ 2 ^ 7 . Experimentally, for Pr < 1 the scaling of 
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Nu versus Ra seems to be independent of Pr, while for Pr < 1 it seems that Nu — Pr a with a small and positive. 
The whole behavior is however not completely clear and it is still under investigation (see |l8| for very recent results). 
At any rate, a negative scaling exponent in the Pr number is ruled out by existing experimental observations. The 
model we have proposed in this paper can partially explain the experimental results. Indeed, the Bolgiano dissipation 
scale rs can be a rather complex function of Pr because either kinematic viscosity or thermal diffusivity can enter in 
the computation of rs- Moreover, if one takes into account intermittent fluctuations and multiscaling effects on the 
dissipation scale, the computation of rg may explain the observed dependence on Pr. This problem deserves more 
investigations in the future. 
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FIGURE CAPTIONS 


Figure 1 Scaling of Nu versus Ra. Data + are taken from |l^| for a 2D numerical simulation (fitted slope 0.280 ± 0.001). 

Data * from our numerical simulations with no-slip boundary conditions (fitted slope 0.283 ± 0.003). Same for 
x but with free-slip boundary conditions on top/bottom walls (fitted slope 0.277 ± 0.003). 

Figure 2 Typical behavior of the energy dissipation (in arbitrary units), e(z ), inside the cell. 

Figure 3 Typical behavior of the Bolgiano length (in lattice units), Lb(z ), inside the cell. 

Figure 4 Scaling of A versus rg. The points corresponds to Ra values of 3.34 • 10 6 ,6.68 • 10 6 ,1.67 • 10'. The straight line 
is a linear fit, yielding A = 1.15 ■ tb- 

Figure 5 Behavior of the fraction of energy dissipated inside the viscous boundary layers over the total energy dissipate 
inside the cell. The line are, respectively a fit with the expected slope 1/7 and with arbitrary slope (~ 0.177). 
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